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Superconductivity of boron-doped diamond, reported recently at T c =4 K, is investigated exploit- 
ing its electronic and vibrational analogies to MgB2. The deformation potential of the hole states 
arising from the C-C bond stretch mode is 60% larger than the corresponding quantity in MgB2 that 
drives its high T c , leading to very large electron-phonon matrix elements. The calculated coupling 
strength A w 0.5 leads to T c in the 5-10 K range and makes phonon coupling the likely mechanism. 
Higher doping should increase T c somewhat, but effects of three dimensionality primarily on the 
density of states keep doped diamond from having a T c closer to that of MgB2 . 

PACS numbers: 71.20.-b, 71.20.Be, 71.20.Eh, 71.27.+a 



Discovery of new types of superconducting mate- 
rials has accelerated since the discovery of the high 
temperature superconductors, with a recent break- 
through being the discovery of superconductivity at 
T c = 40 K in the simple (structurally and electroni- 
cally) compound MgP^.pJ The origin of its remark- 
able superconductivity is now understood to arise 
from charge carriers doped (in this case, self-doped) 
into very strongly bonding states that in turn re- 
spond very sensitively to the bond-stretching vi- 
brational modes. 2] The strong B-B bonds in the 
graphitic layers of MgB2 make it appear near op- 
timal, although graphite itself and diamond are ma- 
terials that have even stronger bonds. Of these two, 
only diamond has bonding states that can conceiv- 
ably become conducting through hole-doping. Q The 
recent report by Ekimov et aL0 of superconduc- 
tivity at 4 K in very heavily boron-doped diamond 
revives the question of mechanisms in strongly co- 
valent materials. Confirmation has been provided 
by Takano et al. who report T c = 7 K in B-doped 
diamond films. 

While the study of B as a hole dopant in dia- 
mond has a long history, there have been recent de- 
velopments due to the ability to dope diamond films 
more heavily (to and beyond the 10 20 cm -3 range) 
than was possible previously. In spite of its growin; 
importance, and unlike the situation for donors 
there has been little theoretical work on the acceptor 
state (such as determining its spatial extent) beyond 
obtaining the structural and vibrational properties 
of the isolated B impurity. An isolated B atom is 
an acceptor with a binding energy of 0.37 eV : 8f for 
which effective mass theory is not applicable, but 
the behavior of B-doped diamond up to and some- 
what beyond the concentration for insulator-metal 
(Mott) transition Cjvf=2xl0 cm -3 has been rather 
well studied experimentally. || (The B concentra- 
tion achieved by Ekimov et al. is about c sc =5xl0 21 
cm = 25cmj with a hole carrier density of nearly 
the same, and introduces a new regime of metal- 



lic diamond that is yet to be understood.) At this 
concentration, B dopants are on average 5-6 A apart 
and their donor states, which form an impurity band 
already at Cm — ^c sc , broaden considerably and 
overlap the valence band edge. In addition, Mamin 
and Inushima have pointed out@ that as the B con- 
centration increases (and they were not yet think- 
ing of the 10 21 range) many of the donor states be- 
come more weakly bound states of B-related com- 
plexes that would encourage formation of broader 
bands. Fontaine has analyzed the concentration de- 
pendence of the activation energyjlfj (0.37 eV at 
low concentration) and concluded that it vanishes 
at c = 8 x 10 20 cm -3 = ^c sc ; for larger concentra- 
tions the system would be a degenerate metal. 

Given these indications of degenerate behavior 
even below 2.5% doping, in this paper we adopt the 
viewpoint that the majority fraction of the hole car- 
riers resides in states overlapping the diamond va- 
lence band, and behave primarily as degenerate va- 
lence band holes. Boeri et aLpH have taken a simi- 
lar viewpoint, and two supercell calculations 
have verified this degenerate metal picture. The dis- 
tinctly different low-concentration, nonmetallic limit 
has also been suggested. We investigate the mag- 
nitude and effect of hole-phonon coupling analo- 
gously to what has been found to drive supercon- 
ductivity in MgB2, and present evidence that at 
hole-doping levels similar to that reported, the hole 
- bond-stretch coupling is surprisingly strong and 
makes phonon exchange a prime candidate for the 
mechanism of pairing. In the case that such cou- 
pling is strong, it can be verified by spectroscopy of 
the Raman-active bond stretch mode. In fact, Eki- 
mov et al. report Q a Raman spectrum in which 
the sharp diamond peak at 1332 cm' 1 has vanished, 
leaving spectral weight in the 1000-1300 cm -1 range. 
This behavior is a plausible extrapolation (consid- 
ering they are very differently prepared materials) 
from more lightly B-doped films in which Ager et al. 
observed an initial weakening and broadening of the 



1332 cm -1 mode,|l5j and a transfer of spectral den- 
sity to peaks in the 940-980 cm' 1 range for concen- 
trations ~ c sc . Zhang et al. reported, for films with 
concentration ~ |c sc a broad peak at 1200 cm -3 
and a very broad feature peaking at 485 cm~ 3 .p^| 

For the sake of definiteness, we consider a hole 
concentration of 0.025/carbon atom, 10% less than 
the B concentration (c sc ) determined for the su- 
perconducting diamond films. At this concentra- 
tion the hole Fermi energy ~Ep lies 0.61 eV below 
the valence band maximum, and the diamond Fermi 
surfaces consist of three zone-centered "spheroids." 
The outer one in particular differs considerably from 
spherical due to the anisotropy of the band mass, 
but effects of anisotropy will be decreased by disor- 
der scattering and in any case would give only second 
order corrections to the properties that we calculate. 

The system we consider is thus 2.5% hole-doped 
diamond. The key points here are (1) the carrier 
states are the very strongly covalent bonding states 
that make diamond so hard, and (2) these states 
should be sensitively coupled to the bond-stretching 
mode, which lies at the very high frequency of 1332 
cm -1 (0.16 eV) in diamond. These ingredients are 
the same as those prevailing in MgB2. There are 
differences, both of a positive and negative nature. 
In MgB2 only two of the nine phonon branches are 
bond-stretching, whereas in diamond these comprise 
three of the six branches. On the other hand, MgB2 
is strongly two dimensional in its important bands (a 
bands), which means a near-step- function increase 
in the density of participating states as doping oc- 
curs; the states in diamond are three-dimensional 
and their Fermi level density of state N(0) increases 
with doping level more slowly. 

A look at the phonon spectrum of diamond [l7| 
reveals that the three optic modes are the bond 
stretching ones, and they have little dispersion so 
0, ~ 0.15 eV is their common unrenormalized fre- 
quency. The theory of carrier-phonon coupling and 
the resulting superconductivity in such systems is 
well developed, and the important features in MgB 2 - 
like systems have been laid out explicitly. The cou- 
pling strength A is given rigorously for an element 
by 

A= E fc ^(Q)<4 2 > _ mjLs m 

M < lo 2 > MujI 1 ' 

where iVf,(0) is the DOS of band b, M is the car- 
bon mass, II =<C \Ib(k,k')\ 2 ~3>fs is the Fermi sur- 
face averaged electron-ion matrix element squared 
for band b, and < u> 2 > is an appropriately defined 
mean square frequency which will simplify to ui 2 , 
the bond stretch frequency renormalized by the hole 
doping. The sum over bands b has been displayed 
explicitly but finally leads to an rms electron-ion ma- 



trix element in the numerator, and the Fermi level 
density of states is N(0)=0.060 states/eV per cell per 
spin. 

Due to non-sphericity and non-parabolicity of 
the three inequivalent bands substantial computa- 
tion would be required to obtain accurate num- 
bers (and anharmonic and non-adiabatic corrections 
would change them, see below). There are two ways 
to obtain approximate values in a pedagogical man- 
ner: (1) calculate the Q=0 deformation potentials 
to obtain the matrix elements for the optic modes, 
or (2) calculate the phonon softening and use the 
lattice dynamical result 

cj 2 q = Q 2 Q + 2Sl Q Re II(Q,0) 
cj 2 = <4^ -> Sl^-2n o A(0) \M\ 2 . (2) 

where n(Q, u>) is the phonon self energy arising from 
the doped holes, and M. is the electron-phonon ma- 
trix element and is determined by I rms (see below). 
We apply both methods to obtain estimates of the 
coupling strength. 

The calculations were done with the Wien2k lin- 
earized augmented plane wave code.[l8j The basis 
size was fixed by R m tK ma x = 7.0 with a sphere ra- 
dius 1.2 for all calculations. While 110 irreducible 
k points were used for pure diamond, 1156 k points 
in the irreducible wedge for the B-doped diamond 
virtual crystal calculations, (nuclear charge Z = 
(1 — x)Zc + xZb = 5.975) because the Fermi sur- 
face volume had to be sampled properly to account 
for screening. Alloy (coherent potential approxima- 
tion, CPA) calculations using the full potential local 
orbital code 19] give bands as in our virtual crystal 
model, the main difference [20j being small disorder 
broadening that would not change our conclusions. 

The central quantity in Eq. ^ is the matrix ele- 
ment, which can be expressed in terms of the de- 
formation potential T>\ we use the definitions of 
Khan and Allen |2l| to avoid ambiguity. V is the 
shift in the hole (valence) band edge with respect to 
the bond stretching motion, whose scale is given by 
u = v / h/2Mfl = 0.034 A. The stretching mode 
is threefold degenerate and can have any direction 
of polarization. We have chosen the polarization in 
which atoms move along a < 111 > direction. Un- 
der this displacement, the threefold eigenvalue splits 
(see Fig. \\\ior the case of doped diamond) at the rate 
of A(e U pp er - e i OW er) k=o / Ad bond = 21 eV/A, where 
dbond is the bond length. Since the twofold band 
splits half as rapidly as the single band (and oppo- 
sitely) this leads to the two deformation potentials 
of magnitude T>\ = 14 eV/A for the nondegener- 
ate band and T>i = 7 cV/A for the doublet, for 
intrinsic diamond. The large deformation potential 
is 60% larger than the (already large) analogous one 
in MgB 2 . 
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FIG. 1: (Color online) Virtual crystal bands of 2.5% B- 
doped diamond without (dashed lines) and with (solid 
lines) a bond stretch phonon frozen in. The atomic dis- 
placement a\^Ax = 0.0309 A is just enough to trans- 
fer all holes to within a single k=0 centered Fermi sur- 
face. The horizontal dashed line indicates the Fermi level 
(aligned for this plot). X denotes the usual zone bound- 
ary point, L designates the zone boundary point in the 
< 111 > direction parallel to the atomic displacements. 



The results for 2.5% B-doping are needed for cal- 
culation of the coupling strength, and are shown 
in Fig. They are renormalized by B-doping by 
3% (downward) from those of intrinsic diamond, so 
again we have V x =14 eV/A, V 2 = 7.0 eV/A. These 
deformation potentials are undoubtedly the largest 
yet encountered for any metallic solid, being directly 
related to the great bond strength of diamond. Since 
the three deformation potentials contribute addi- 
tively to the coupling strength, we simplify by using 
the root mean square value I rm s = 10 eV/A. The 
rms electron-phonon matrix element, to be used be- 
low, is M. = \/w o /fi u I rma = 0.70 eV; here u is 
the renormalized optic frequency. 

Together with the value N(0) = 0.060 states/eV- 
spin, MVL 2 = 65 eV/A 2 , and (calculated below) 
uo 2 = O.68f2o, we obtain from Eq. fl] the coupling 
strength A = 0.55. The coupling is confined to a 
set of three optic branches which comprise a narrow 
peak centered at u . The conventional theory, ne- 
glecting very minor strong-coupling corrections. |22j 
gives T c = (wo/1.2) exp[-l/(j^ - (J,*)], where fi* 
is an effective Coulomb repulsion that is uncertain 
for doped diamond. Using the conventional value fi* 
= 0.15 with uj = 0.128 eV gives T c = 9 K, gratify- 
ingly (and probably fortuitously) close to the exper- 
imental values of 4 K to 7 K. To obtain the initially 
observed value T c =4 K would require A= 0.48, or al- 
ternatively n* « 0.20, i.e. relatively small changes. 

A less direct way of obtaining the coupling 
strength, but one that (numerically) includes averag- 
ing over bands properly, is to calculate the renormal- 
ized phonon frequency and apply Eq. |2 The calcu- 
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FIG. 2: (Color online) Plots of the energy of dis- 
tortion for the frozen-in bond stretch mode, for (top 
straight line) undoped diamond, and (bottom line) 2.5% 
B-doped diamond. The coordinate uc is for one of 
the two identically displaced atoms. Pristine diamond 
follows a quadratic plus lowest-order anharmonic form 
AE(u) — A2U 2 + A3U 3 accurately, as indicated by the 
straight dashed line. After doping the AE(u) functional 
form becomes very complex. The horizontal arrows in- 
dicate the atomic displacements for which one (or two) 
Fermi surfaces disappear. The inset shows the AE(u) 
curves themselves. 



lated change in energy versus atomic displacement, 
plotted as AE(u)/u 2 , is shown in Fig. |2 both for 
intrinsic and doped diamond. The difference due to 
doping is striking. The result for diamond is simple 
to understand: the harmonic u 2 term give s Qharm 
= 1308 cm -1 similar to literature values. [171 l23l] and 
the A3U 3 term quantities its anharmonicity. 

The AE(u)/u 2 curve for the doped case (see Fig. 
EJ) is much more complex. The reason is clarified by 
the red arrows on the plot, which mark the displace- 
ments where some piece(s) of Fermi surface vanishes. 
One of these values of displacement is also used for 
the deformation potential plot of Fig. ffl where one 
shifted band edge is lying exactly at Ejr. Boeri et al. 
have described how, at such topological transitions 
of the Fermi surface, the energy is non-analytic. [24j 
In Fig. |21 one can imagine a straight line behavior 
similar to (and nearly parallel to) that of diamond 
for uc between the two topological transitions, with 
changes of behavior occurring beyond each transi- 
tion point. It is of special relevance that these po- 
sitions are roughly at the bond-stretch amplitude, 
hence they are physically important. 

Returning to the coefficient of the u 2 (harmonic) 
term, for the doped case it is 0.68 that of diamond, 
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giving (in harmonic approximation) uj = 1070 cm . 
This renormalization of the square of the optic mode 
frequency (see Eq. [2Jl by 32% is a vivid indication of 
the strong coupling, even for the case of 2.5% holes. 
Substituting ui ratio into Eq. fallows us to extract 
the electron-phonon matrix element Ai = 0.67 eV, 
within 5% of the value obtained from the deforma- 
tion calculation. Combined with the deformation 
potential result, the predicted coupling strength is 
A = 0.53 ± 0.03. As mentioned above, this value is 
quite consistent with the observed critical tempera- 
ture, but certainly such good agreement may not be 
warranted. 

Our treatment neglects some complicating fea- 
tures. The Jahn- Teller splitting of the isolated 
B substitutional impurity (0.8 cm -1 from Fabry- 
Perot spectroscopy [25| ) is three orders of magni- 
tude smaller than energy differences involved in the 
bond-stretch mode and therefore is negligible. It 
has been suggested^, El that supercells offer a 
more realistic model than a virtual crystal treat- 
ment. Our calculation of a C31B supercell indi- 
cated strong ordered-boron effects (which are un- 
physical), and our CPA alloy calculations!^ give 
bands like the virtual crystal model, with small dis- 
order broadening added. Another factor is anhar- 
monicity, which includes conventional anharmonic- 
ity and the non-adiabatic effects that cause the non- 
linearity of AE(u)/u 2 in Fig. |2 Making the an- 
harmonic corrections need not change the effective 
phonon frequency greatly, as shown for MgB2 by 
Lazzeri et aL[2(| who found that for MgB 2 three- 
and four-phonon corrections gave strongly canceling 
corrections to the vibrational frequency. The valid- 



ity of the Migdal-Eliashberg itself becomes an inter- 
esting question, and more so for lower doping levels. 
For the 2.5% concentration considered here, the ra- 
tio of phonon frequency to electron energy scales is 
u>/Ep = 0.25, certainly not the small parameter that 
is usually envisioned as a perturbation expansion pa- 
rameter. Doped diamond provides a new system in 
which to investigate non-adiabatic effects. 

Now we summarize. Based on the experimen- 
tal information available so far, the B doping level 
in diamond achieved by Ekimov et al. should re- 
sult in hole doping of the diamond valence bands 
to a level Ej? w 0.6 eV. Calculations bear out the 
analogy to MgB2 that deformation potentials due to 
bond stretching are extremely large, and evaluation 
of the holc-phonon coupling strength using conven- 
tional theory leads to A w 0.55, a renormalization 
of the optic mode frequency by -20%, and T c in the 
5-10 K range. These results indicate that phonon 
coupling is the likely candidate for the pairing mech- 
anism, consistent with the conclusions of Boeri et 
q/.pr| The low carrier density (for a metal) implies 
both poor screening of the Coulomb interaction and 
the intrusion of non-adiabatic effects, which are pri- 
mary candidates for further study. Higher doping 
levels should increase T c , but probably not to any- 
thing like that occurring in MgB2. Definitive calcu- 
lations will require CPA calculations of the electron- 
phonon coupling characteristics. 
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